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Abstract 

As is well known from the work of R. Glassey and J. Schaef- 
. fer the main energy estimates which are used in global existence 

^ | results for the gravitational Vlasov-Poisson system do not apply to 

$h ' the relativistic version of this system, and smooth solutions to the 

initial value problem with spherically symmetric initial data of nega- 
tive energy blow up in finite time. For similar reasons the variational 
techniques by which Y. Guo and G. Rein obtained nonlinear sta- 
bility results for the Vlasov-Poisson system El El E3 GB1 EES] do not 
apply in the relativistic situation. In the present paper a direct, non- 
variational approach is used to prove nonlinear stability of certain 
steady states of the relativistic Vlasov-Poisson system against spheri- 
cally symmetric, dynamically accessible perturbations. The resulting 
stability estimates imply that smooth solutions with spherically sym- 
metric initial data which are sufficiently close to the stable steady 
states exist globally in time. 
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1 Introduction 



The topic of the present investigation is the following nonlinear system of par- 
tial differential equations, known as the relativistic Vlasov-Poisson system: 

d t f + ^J==.V x f-V x U-V v f = 0, (1.1) 

v 1 + M 

AU = 4np, lim U(t,x) = 0, (1.2) 

|a;|— >oo 



p(t,x)= J f(t,x,v)dv. (1.3) 

Here teM denotes time, x,t> GlR 3 denote position and momentum, and 
f = f(t,x,v) >0 is the time-dependent density on phase space of a large en- 
semble of particles which interact only by the Newtonian gravitational po- 
tential U = U(t,x) which the ensemble creates collectively through its spatial 
mass density p = p(t,x). Collisions among the particles are assumed to be 
sufficiently rare to be neglected. All particles are assumed to be of the same 
rest mass, and all physical constants such as the rest mass of a particle, the 
gravitational constant, and the speed of light are normalized to unity. The 
individual particles obey the following equations of motion, which form the 
characteristic system of the Vlasov equation (jl.ljl : 

v 

x= f=== , v = -V x U(t,x). 



The system is called relativistic because of the relation between the momen- 
tum v of a particle and its velocity x. However, the system as a whole is 
neither Lorentz nor Galilei invariant. 

In astrophysics galaxies or globular clusters are often modeled as large en- 
sembles of particles, i.e., stars, in this way. But from the point of view of such 



applications the non-relativistic Vlasov-Poisson system where v/y/l + \v | 2 is 
replaced by v in the Vlasov equation is much more important, while a truly 
relativistic formulation is given by the Vlasov-Einstein system PQ. From a 
mathematics point of view the non-relativistic Vlasov-Poisson system is by 
now quite well understood. In particular, the initial value problem has global 
weak and, for smooth initial data, global classical solutions, cf. [TU [T71 ■ 
Moreover, nonlinearly stable steady states have been established by varia- 
tional techniques, i.e., by minimizing suitably chosen energy- Casimir func- 
tionals, cf . [UJ |H1 El Ed UH1 EES] • The stability question has also received 
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a lot of attention in the astrophysics literature; we refer to [3J 0] and the 
references there. An extensive review of the mathematical results on global 
existence and stability for the non-relativistic system can be found in [2*U] . 

For the relativistic version of the system stated above both the global 
existence results and the variational approach to stability fail, and we briefly 
explain why The total energy 

E u M + E pot (f)--= J J V± + Wf(x,v)dvdx-^ J \VU f (x)\ 2 dx 

of a state / with induced potential Uf is conserved along solutions, and 
so is the mass ||/||i or any other L p norm ||/|| p . Whether one wishes to 
minimize the energy in a variational stability analysis or to extract an a- 
priori bound from the energy for a global existence result, in both cases 
one needs to control the negative potential energy in terms of a power less 
than one of the kinetic energy and the other conserved quantities. But the 
Hardy-Littlewood-Sobolev inequality together with standard interpolation 
arguments imply that 

-^ P ot(/)<c'ii/iir fe)/3 ii/ii^ /3 J E km (/), 

where < k < 2, so no matter which k we choose, in the context of the above 
estimate E pot (f) is of the same order as E kin (f). This is the reason why 
solutions can blow up in finite time, as shown in j3] for spherically symmetric 
initial data of negative energy, and why the variational stability approach 
fails. For the non-relativistic system \v | 2 /2 replaces -y/l + \v \ 2 in the kinetic 
energy, and as a consequence the latter appears with the exponent 1/2 in the 
above estimate. 

Recently a non-variational stability approach was introduced for the 
Vlasov-Poisson system ^U|- Its advantage is that it can deal with steady 
states which are only local minimizers of an energy- Casimir functional and 
not global ones, in particular, it requires no lower bound on the energy- 
Casimir functional. For the reason explained above this becomes essential in 
the case of the relativistic Vlasov-Poisson system, and it is the aim of the 
present investigation to show that this new method applies to a variety of 
steady states of the relativistic Vlasov-Poisson system and provides their sta- 
bility against spherically symmetric, dynamically accessible perturbations — 
to remove the symmetry assumption is an open problem. In addition we show 
that the resulting stability estimates provide bounds on spherically symmet- 
ric solutions starting close to a stable steady state, which are sufficient to 
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prove global existence. To our knowledge the resulting global existence re- 
sults are the first such results for initial data which are not small. 

The paper proceeds as follows. In the next section we state and explain 
our main results. Section 3 contains the stability analysis. In the Newtonian 
case it has been known for a long time in the astrophysics literature that 
the quadratic term in an expansion of a suitable energy- Casimir functional 
is indeed positive definite on the set of linearly dynamically accessible per- 
turbations. We prove such a result, which can be used to study linearized 
stability, for the relativistic Vlasov-Poisson system. The essential difficulty 
is then to connect this linear result to the nonlinear problem. Here our anal- 
ysis proceeds essentially as in [TU], but we strive for greater generality as 
far as the admissible steady states are concerned, because as opposed to the 
Newtonian case no other stability results are available here. In Section 4 we 
briefly show how the stability estimates imply global existence for the corre- 
sponding, perturbed solutions. In the last section we present a large variety 
of steady states to which our method applies. So far only the existence of 
the so-called polytropic steady states has been established for the relativistic 
Vlasov-Poisson system j2]. 

2 Main results 

Consider a steady state fo of the relativistic Vlasov-Poisson system with 
induced potential Uq and spatial density p , where f is a function of the 
particle energy: 

f (x,v) = 4>(E) where E = E(x,v) := y/l + \v\ 2 + U (x). (2.1) 

For a time-independent potential Uq the particle energy E is conserved along 
characteristics of the Vlasov equation so that any function of the particle 
energy satisfies the Vlasov equation with potential Uq, and it remains to 
show that for a given choice of the Poisson equation for Uq has a solution, 
where the spatial density now becomes a functional of Uq. In Section |3] we 
present a variety of functions for which this approach gives a compactly 
supported steady state with finite mass. Steady states of the form (|2.1j) 
are called isotropic, and it can be shown that they are always spherically 
symmetric, i.e., 

f (x,v) = f {Ax,Av), AgSO(3), x,veR 3 . 



4 



We want to analyze the stability of such isotropic steady states. The total 
energy 

H(f):=E km (f)+E pot (f) = JJ ^l + \v\ 2 f(x,v)dvdx-±- J \VU f (x)\ 2 dx 

of a state / is conserved along smooth solutions of the relativistic Vlasov- 
Poisson system. But if we expand TC about any state / with potential U 
we find that 

W) = Wo) + JJ (VT+W+Uo) (f-f )dvdx-^J \VU f -VU \ 2 dx, 

and the linear part in the expansion does not vanish. Hence we cannot 
use the energy as a Lyapunov function in a stability analysis. To remedy 
this situation we observe that for any reasonable function $ the Casimir 
functional 

C(/):= [[ $(f{x,v))dvdx 



is conserved as well. If we expand the energy-Casimir functional 

n c -.=n+c 

about an isotropic steady state, then with E defined as in Eqn. (|2.1|) . 

naif) = n c (fo) + JJ( E + $ U)) (/ - fo) dvdx 

/ \VU f -VU \ 2 dx+ l - JJ$"(f )(f-f ) 2 dvdx + .... 

We want to choose <£> is such a way that at least formally fo is a critical point 
of the energy-Casimir functional, i.e., Q'(fo) = — E. In order to make this 
rigorous we specify the following assumptions on <ft and /o : 

Assumptions on <p and fa: (a) 0eC7(R), there exists a cut-off energy Eq 
such that <f)(E) = for E>E , lim £ _„ oo 0(E) = oo, and G C 2 (] - oo,E [) 
with 

< for E < E Q , and liminf cj)'(E) > -oo. 



E—*Et 



o- 



(b) /o£Ce(U 6 ) is compactly supported, satisfies the relation (|2.1|) . where U 
denotes the potential induced by f , and is non-trivial. 
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Since 

supp/o = {(x,v) G R 6 | E(x,v) = y/i + \ v \2 + U (x) < E } 

and since Uq is spherically symmetric and strictly increasing as a function 
of r = |x| with lim r ^ 00 ?7o( r ) = 0, part (b) necessarily implies that E <1. In 
the last section we provide a large class of steady states which satisfy these 
assumptions. Examples are polytropic steady states 

<i>(E) = (E -E) k + , 

with 1 < k < 7/2, and the King model 

<P(E) = {e E °- E -l) + ; 

(•)+ denotes the positive part. 

For a function satisfying the assumption above, (p:] — oo,E ] 
is continuous and invertible with continuous inverse _1 : [0,oo[— >] - 
and we define $ : [0, oo[^ R by 

$(/):=- f <j>-\z)dz, /g[0,oo[. (2.2) 
Jo 

In particular, $ G C 1 ([0,oo[). 

A crucial step in the stability analysis is to specify the set of admissible 
perturbations. From a physics point of view perturbations arise by some 
exterior force field acting on the ensemble represented by the steady state 
/o- Such a field induces a measure preserving flow on phase space which 
redistributes the particles. We refer to perturbations of the form f — fo° 
T with T : M 6 — > M 6 a measure preserving C 1 -diffeomorphism as dynamically 
accessible from f . For technical reasons we have to restrict ourselves to 
spherically symmetric such perturbations. More precisely, we say that the 
C 1 -diffeomorphism T : M 6 — > M 6 respects spherical symmetry if for all x,v G 1R 3 
and all rotations A<E SO (3), 

T(Ax,Av) = (Ax',Av') and \x' x v'\ — \x x v\, where (x',v')—T(x,v). 

Such a redistribution of particles on phase space would be caused by the 
action of a field which is spherically symmetric, and from a physics point of 
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view this restriction is undesirable. The set of admissible perturbations is 
defined as 

Vf Q :— j/ = / oT | T:M 6 — >IR 6 is a measure preserving C 1 -diffeomorphism 

which respects spherical symmetry |. 

It is important to note that this set is invariant under classical solutions of 
the relativistic Vlasov-Poisson system. We are going to measure the distance 
of a state f &Vf from the steady state fo by the quantity 

d(/,/o):= JjMf)-Hf )+E(f-f )]dvdx+±- J \VU f -VU \ 2 dx, 

which is closely related to the energy-Casimir functional: 

d{f,f )=Hc(f)-Hc(fo) + ^ J \VU f -VU \ 2 dx. (2.3) 

As we will see in the next section there exists a constant C > which depends 
only on the steady state fo such that 

1 1/ - /ol II + 1 1 VI// - VI7 | |1 < C7d(/, /o), feVfy. (2.4) 

The major part of the analysis will be concerned with proving the following 
result which says — in a precise, quantified manner — that the steady state is 
a local minimizer of the energy-Casimir functional in the set T>f . 

Theorem 2.1 There exist constants 5o>0 and Co>0 such that for all f e 
Vf with d(f,fo) <5o the following estimate holds: 

H c (f)-Hc(f )>C \\VU f -VU \\l 

This result is proven in the next section. Our main result is the following 
theorem, which is an immediate corollary. 

Theorem 2.2 There exist constants 5>0 and C > such that for any initial 

o 

datum f&Vf with 

d(f,fo)<5 

the corresponding solution ti—>f(t) of the relativistic Vlasov-Poisson system 

o 

with /(0) = / exists globally in time and satisfies the estimate 

d(f(t)Jo)<Cd(fJ ), t>0. 
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Proof. Let 5 := 5 (1 + l/(47rC )) _1 with S and C from Theorem O Con- 
sider a solution [0,T[3t*— > f(t) of the relativistic Vlasov-Poisson system with 

o 

fET>f on some maximal interval of existence; as to such a local existence 
result we refer to the comments in Section Now assume that 

d(f,f )<5<5 . 

By continuity we can choose some maximal t* G]0,T] such that 

d(f(t),f )<5 , te[o,t% 

Now f(t) eDj for all t G [0,T[, and hence Theorem 12.11 the relation (|2.3|) of 
<i to the energy-Casimir functional, and the fact that the latter is a conserved 
quantity yield the following chain of estimates for tE [0,t*[: 

d{f(t),f ) = H c {f{t)) - Hc{fo) + ^ 1 1 VU m - VC/o| \l 

< Hc{f{t))-H c {h) + -^r (H c (f(t))-H c (fo)) 

= { 1 + ^C-) ( ^ (/(0)) - HM)) * {} + 4^) d ^ /o) < *>■ 

This implies that t*=T. By Proposition 14. 1[ T = oo, and Theorem 12.21 is 
established. □ 

Remark. Theorems 12. II and 12.21 are analogous to the corresponding results 
obtained for the non- relativistic case in ^U], but there are two differences. 
Firstly, the analysis in |10j is restricted to the King model, the reason being 
that stability for other models had already been obtained via variational 
methods. Since these methods do not work in the relativistic context, as 
was explained in the introduction, we keep the admissible steady states more 
general here, but the same generality would have been possible in the non- 
relativistic case. Secondly, for the non-relativistic case global existence of 
classical solutions is known for general data, independently of any stability 
analysis. In the relativistic case this is not so. Theorem 12.21 includes a 
new global existence result for the relativistic system for data which are not 
subject to a size restriction. This should be compared with the results in 
jS], in particular, with the fact that spherically symmetric solutions with 

o 

TC(f) <0 blow up in finite time. In particular, Theorem 12.21 answers an open 
problem which was formulated in j2J Final remarks]. 
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3 Proof of Theorem EH 



The essence of the proof of Theorem 12. II is the analysis of the quadratic term 

D 2 H G {h)[g]:= l ~ [[ &'(f )g 2 dvdx-±- [ \VU g \ 2 dx (3.1) 
2 JJ{f 0>0 } 8tt J 

which arises in the expansion of the energy- Casimir functional Tic- If Theo- 
rem were false, then a tangent direction g to the set T>/ of dynamically 
accessible states would exist on which this quadratic part is negative, cf. 
Lemma 13.21 On the other hand, if we define the Poisson bracket of two 
functions f,h:R 6 ^R by 

{f,h}:=VJ-V v h- V v f ■ V x h, (3.2) 

then for linearized, dynamically accessible states g = {fo,h} one can show that 
the quadratic term is indeed positive definite, cf. Lemma l3~3| we do not go 
into the symplectic dynamics details behind this terminology and construc- 
tion. Note however that at least formally, states of this bracket form arise as 
tangent vectors to the manifold T>f at the point fo, and that the set of these 
states is invariant under the linearized system. The analogue of Lemma 13.31 
for the non-relativistic case is well known in the astrophysics literature, where 
it has been used to prove linearized stability [121 HU I2I| , but for the relativis- 
ts case the result seems to be new. To conclude the proof of Theorem 12 .11 we 
then show that the function g with D 2 7ic(fo)[g\ <0 obtained in Lemma f3. 21 
can be written in the form g = {fo,h}, which contradicts Lemma 13.31 At 
some points the arguments are similar to the non-relativistic case considered 
in |1U| . But since we consider a different system and also a general class of 
admissible steady states and not just one example, the analysis is technically 
more difficult, and we prefer to give a self-contained proof for the present 
situation. 

At several points we have to exploit the spherical symmetry assumption 
which is part of our definition of the set T>f of dynamically accessible data. 
To do so we on occasion use coordinates which are adapted to the spherical 
symmetry: 

r:=|x|, w:= , L:=|xxt>| 2 ; (3.3) 

r 

w is the radial component of momentum and L is the modulus of angular 
momentum squared. It is easy to see that for any spherically symmetric 
function / by abuse of notation, f(x,v) = f(r,w,L). 
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We first collect some properties of the function $ defined in Eqn. (J2.2I) . 

Lemma 3.1 (a) $ G C^^ooQ nC 3 (]0,oo[) ; $(f)>-E f for f>0, and 
for / > 0, 

'I''!/) o '(/). ^(/) = - — J— - , $"'(/)- -" ( ' rJ( ^ ) 



(b) T/iere exists C > swc/i t/iai /or a// functions f <EVf , 

JJ mf)-$(f ) + E(f-f )}dvdx>C JJ \f-f \ 2 dvdx, 

in particular, Eqn. \2.4\j holds. 

Proof. The formulas for the derivatives are obvious, and the lower bound 
for $ follows from the fact that < E . The assumptions on <// imply that 

C:=^inf{$"(/)|0</<||/o||oo + l}>0. 

For feV f0 , 

E(f - / ) = -&(f Q ) (/ - /o) on {/ > 0} 

while 

E(f-f o )=Ef>E o f = -<S> f (0)f = -&(f o )(f-f o ) on {/ = 0}. 
Hence by Taylor's Theorem 

*(/)-*(/o)+£(/-/o) > ^(/)-$(/o)-^(/o)(/-/o) 

= lim [$(/ + e ) -$(/ + e )-^(/ + e )(/-/ )] 

= hm J$"fe)|/-/o| 2 >^|/-/o| 2 , 

e— >0+ Z 

since < £ e < 1 1/ | |oo + 1 ase^O, and the proof is complete. □ 
Assuming that Theorem 12.11 were false we construct a state g on which 
the quadratic term D 2 TCc{fo) is negative, more precisely: 
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Lemma 3.2 Assume that Theorem \2. 1\ were false. Then there exists a func- 
tion g£L 2 (M. 6 ) which is spherically symmetric, supported in supp/o? even in 
v, i.e., g(x,—v)=g(x,v), and such that 

^l|W ff ||^ = l, (3.4) 



&Hc{hM = \ II <5>"(f )g 2 dvdx-l<0, (3.5) 

Z JJ{fo>0} 

and for all functions G = G(f,L) G C 2 ([0,oo[ 2 ) with G(0,L) = d f G(0,L) = 
for L>0 and djG bounded, 

J J d f G(f ,L)gdvdx = 0. (3.6) 

Remark. Eqn. (|3.6|) expresses the fact that g is a direction tangent to Vf 
at f . If ti— >f(r) is a curve on Vf with f(0) = f and 4-f(0)—g then by 
definition of Vf , JJ G(f(r),L) = JJ G(fo,L) so that the formal derivative of 
the left hand side at r = vanishes, which is Eqn. ()3.6|) . 
Proof of Lemma 13. 2L Since we assume that Theorem 12.11 is false, there 
exists a sequence (f n ) cP/ such that for all neN, 

d{fn,fo)<-, 

n 



but 



If we let 



nc(fn)-H (fo)<^-\\VU fn -VU \\l (3.7) 



57m 
1 



i.e.. 



fn = fo + cr n g n with — ||VC/ 9n ||2 = l, (3.8) 

1 1 

^n:=^=||W/„- Wo||2, 9n-=—(fn-fo), 

V87T cr n 

then in particular, 

a 2 n <d(f n ,f )<-. (3.9) 
n 
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A weak limit g of a subsequence of (g n ). 

Using the definition of d and Eqns. (|3.7j) . (|3.8|) . and (|2.3j) we find that 



1 



[<$>(f n )-$(fo)+Ea n g n ]dvdx-l 
= \(d{f n ,fo)-h\VU fn -VU \\l 

= l(7fc(/.)-« (/.))<~llv^-vn ) |B=i. 

By Lemma f3. II this implies that 

l + ^>^ JJ mf n )-$(f )+Ea n g n }dvdx>C J J \g n \ 2 dvdx, (3.10) 

and hence the sequence (g n ) is bounded in L 2 (M 6 ). We extract a subsequence, 
again denoted by (g n ), such that 

gn^g weakly in L 2 (R 6 ). 

Since the functions g n are spherically symmetric so is g, and we need to show 
that it is supported on supp/ - The term in brackets in Eqn. (|3.1U|) is non- 
negative, and on the set {E > E } the steady state distribution f and hence 
also $(/ ) vanish, while ^{fn)>-E fn = -E o- n g n . Hence 

2>^/7 Mfn)-$(fo) + Ea n g n ]dvdx>— [[ (E — E ) g n dvdx, 

°~nJJ{E>E } a nJJ{E>E } 

and by 



// 

J J{. 



(E — E ) g n dvdx <2a n — >0, n— >oo. (3-H) 

{E>E } 

For any fixed Eq < E\ < 1 , 

ff E~E 

g n dvdx< / / — — g n dvdx — >0; (3-12) 

{E>E 1 } JJ{E>E } El — Eq 

notice that g n > on the set {E > Eq}. Since ()3.12|) holds for any Eq < E\ < 1, 
we conclude that g is supported on supp/o = {E< E Q }. Moreover, since for 
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E <Ei<l fixed the set {E<E{\ is compact, the L 2 -bound on (g n ) and 
Eqn. (J3.12j) imply that (g n ) is bounded in L 1 (IR 6 ) as well. 

Proof of jg. 

First we note that Uq is spherically symmetric and radially increasing, so in 
particular, Uq(0) <U <0. Hence 



1 + \v \ 2 \g n \ dv dx < // E\g n \dvdx — Uq(0) // \g n \dvdx 



< E \g n \dvdx+ // (E — E )\g n \dvdx 

J J JJ{E>E } 

-U (0) / / \g n \dvdx, 



so that together with (j3.11|) and the L 1 -bound on (g n ) the kinetic energy 
is bounded along (g n ). By well known interpolation arguments |2"Ul Ch. 1, 
Lemma 5.1] the sequence of induced spatial densities (p 9n ) is bounded in 
L 5//4 (M 3 ), so without loss of generality this sequence converges weakly in 
L 5 / 4 (M 3 ). We fix E <E 1 <1 and such that 1 + U (R 1 ) = E 1 , which 

implies that E(x,v) >E\ for \x\ >R\. Then by (J3.12j) . 



\p g Jdx< g n dvdx—>0, 
li} JJiEyE!} 



'{\x\>Ri 

i.e., the sequence (p 9n ) remains concentrated, and hence 

VU gn ^VU g strongly in L 2 (M 3 ), 

cf. pUl Ch. 2, Lemma 3.2]. Passing to the limit in Eqn. ()3.8|) proves that g 
satisfies the Eqn. ()3.4I) . 

Proof of UHXl . 

We first claim that there exists a sequence of sets Kj C Kj + i C ... C supp fo 
such that for a subsequence which we again denote by (g n ), 

vol (supp/o \Kj) <- and lim a n g n = uniformly on Kj, j G N. 

j n^oo 

To see this we note that | |cr n p n | I2 < Ca n — >• so that a subsequence converges 
to zero pointwise a. e.. Then we use Egorov's theorem on the set supp/o 
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which has finite measure. In the following arguments we need to stay away 
from the boundary of the latter set, and hence we define 

S m := { (x, v) G M 6 | E{x, v) < E - 1/m} , m G N. 

Clearly, 5 m := inf s m fo > 0. Hence for all sufficiently large n, 

S m /2<f + a n g n < U/olloo + 1 on S m nKj. 

On the set S m C\Kj, 

$(/n)-$(/o)+Wn-/o) = ^'U)^^ 

A O 

with 0<£< 1; recall that $ G C 3 (]0,oo[). In particular, 

\^"(f + ^ n g n )\<sup{\^"(z)\ | 5 m /2<z< ||/ ||oo + l}=:C m <oo. 



Using (|3.10|) we find that for all sufficiently large 



n, 



i U &'{f Q )\g n \ 2 dvdx=± [[ Mf n )-*{fo)+E(f n -f )]dvdx 

Z J JS m r\Kj °n J J S m DKj 

~£~2 &"(fo+frn9n){Vn9n) 3 dvdx 
bcr n JJSmnKj 

< 1 + - + C m sup|cr„5f n | // \g n \ 2 dvdx. 
n k 3 JJ 

Now g n ^g weakly in L 2 (M 6 ) and o n g n — >0 uniformly on Taking the 
limit n — > oo implies that for all j,m G N, 



z JJs m nK4 



2 

the latter expression is L 2 -weakly lower semicontinuous in g. By the bound- 
edness of <&"(/o) on and the choice of Kj we find with j — > oo that 



i/7 $"(/ y^^<i, 



2 

and with m— >oo the monotone convergence theorem implies Eqn. (|3.5|) . 
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Proof of arm . 

Let G = G(f,L) be a function as specified in Eqn. (|3.6j) . By Taylor expansion 
with respect to the first argument, 

G(f n ,L) - G(f ,L) = d f G(f ,L) a n g n + -d 2 f G(f + ra n g n ,L) {a n g n ) 2 

for some rG [0,1]. For f eV fo , jj G(f,L) = jj G(f ,L), and hence 



d f G(f ,L)g n dvdx = --o- n // d 2 f G(f + Ta n g n ,L) g 2 n dvdx ^0; 



'f'^\.JO,^)gnU'VUX = --U n J J U f ( 

note that d 2 G is bounded, (g n ) is bounded in L 2 (M 6 ), and a n — >0. On the 
other hand dfG(fo,L) is supported on the compact set supp/o and hence 
bounded. Since g n ^9 weakly in L 2 (M 6 ), Eqn. (|3.6|) follows as n^-oo. 
Conclusion of the proof of Lemma \3. 6 A 

The function g constructed above has all the required properties, except that 
it need not be even in v. However, if we decompose it into its even and odd 
parts with respect to v, g = g eve n + SWd, then g even satisfies (|3.4|) . (|3.5|) . (|3.6|) 
as well. Since p g = p gcvcn we have V£/ Seven = VU g , and ()3.4|) remains valid. 
Since $"(/o) > is even in v, 



2 



1>U[ <&"(/o) ^even + Soddf = \ ff $"(/ ) ((^ven) 2 + (Sodd) 2 ) 

>lf[ $'U) (Seven) 2 , 

'{/o>0} 



i.e., (|3.5|) remains valid. Finally, for G as in (|3.6|) . dfG(f ,L) is even in v so 
that the odd part of g drops out of Eqn. (|3.6|) . and the proof of Lemma [3.21 
is complete. □ 
As was said above for states of the form g = {fo,h} the quadratic term in 
the expansion of Tic is positive definite; the Poisson bracket was defined in 
Eqn. ()3.2j) . It turns out that it will make some arguments technically easier 
later on, if we prove this fact for states of the form g = {E,h}. 

Lemma 3.3 Let h G C^°(M 6 ) be spherically symmetric with supp/i C {/o > 0} 
and such that h(x,—v) = —h(x,v). Then the following inequality holds: 

D 2 Hc(fo)[{E,h}} 

1 ff 1 



> -- 



2JJ <j>'(E) 



\x-v\ 



E,— 



x-v 



\. U'h 2 



r(l + |tf) 3 / 2 



dvdx. 
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Proof. Let 



%(*):=// 



{EM 
\x-y\ 



dvdy 



denote the potential induced by —{E,h}. Using the definition of the Poisson 
bracket, 

/{E,h}dv = 'V x - / =h(x,v)dv. 
J a/i + I^I 2 

Since both h and Uh are spherically symmetric, 

w 



U' h (r)=4n 



71" 



V 



h(x,v)dv. 



By the Cauchy-Schwarz inequality, 



1 

8^ 
Since 



\VU h \ 2 dx<2ir 



-w 



b\E)dv 



w 



•\E)=w^<j ) ( KV /l+w 2 + L/r 2 + U (r^=w-^ ( j ) (E), 



d 



dv 



dx. 



an integration by parts with respect to w yields 



w 



^E)dv = p (r):=p fQ (r) 



Hence 



D^ c (h)[ { EM\>-\fj^ 



(E) 



v 1 + |f I 



dvdx. 



Since h is odd in v and thus in w the function 

u,(r,w,L) := — h(r,w,L) 
rw 

is smooth away from r — 0. Using the identity h = rwp a straight forward 
computation shows that 



\{E,h}\ 



2 / \2 

= (rw) 



\{E,fj}\ 2 + {E^ 2 rw{E,rw}}-p 2 rw{E, {E,rw}}. 
16 



The first term is as claimed in the lemma. The second term leads to 
{ip{E) , fi 2 rw{E ,rw}} with ip a primitive of 1/0', and the integral of this 
expression with respect to x and v vanishes via an integration by parts; if we 
cut a small ball of radius e about x = from the x-integral then the surface 
integral appearing after the integration by parts with respect to x vanishes 
for e— >0 since rp?<C/r. Finally, another lengthy, but straight forward 
computation using the fact that AUo = Anp implies that 

2 r „ rr, ^ 1 h 2 m(r) , x h 2 

-fi 2 rw{E,{E,rw}}= - ■ , ^l + ^p^r) 



1 + M 2 v/l + M 2 r 'y/l + W' 

If we substitute this into the above estimate for D 2 7ic the proof of Lemma EOl 
is complete. □ 

In order to derive a contradiction between Lemma 13.31 and Lemma 13.21 
we have to show that the function g provided by the latter can be written in 
the form {E,h}. Taking the opposite sign is of course equivalent but more 
convenient in the sequel. 
Solving g = {—E,h}. 

Let g be the function obtained in Lemma 13.21 In the spherical variables 
introduced in (j3.3J) the equation to be solved for h takes the form 

d r h+[ = -%{r) \d w h = g. (3.13) 



\Jl-\-w 2 + L/r 2 \r 3 ^l + w 2 + L/r 
The characteristic system of this equation reads 



w L . , . , 

**= 5 7, . , . r r^-Uo(r). (3.14) 



v/l + w 2 + L/r 2 ' r i ^/l + w 2 + L/r 2 

In order to analyze it, we introduce for fixed L> the function 



= U (r) + y/l + L/r 2 . 
In the spherical variables the particle energy takes the form 



E = E(x,v)=E(r,w,L) := y^l + w 2 + L/r 2 + U (r). (3.15) 

It is conserved along solutions of the characteristic system IJ3.14J) . while L 
only plays the role of a parameter. For given values of E and L we want to 
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identify the r-interval in which a corresponding characteristic can range. By 
spherical symmetry, 

^o( r ) = 2 — ' wnere m o( r ) : = 47r / p (s)s 2 ds, r>0. 

r Jo 

Since the steady state is non-trivial and U increasing, U (0) < E — 1, and 
by (|2DJ, Po(0)>0. Next, 

*'lW=0 ^ m o (r)-^=£= = 0, 

y L + r z 

and since the left hand side of the latter equation is strictly increasing for 
L > with limit — \f~L for r — > and | |/o| |i > for r — > oo, there exists a unique 
radius rx > such that 

^' L (r L ) = 0, ^(r) < for r < r L , ^' L (r) > for r > r L . 

Since lim r _ >0 x I / L( r ) = oo and lim r ^ 00 \l/i(r) = 1 we conclude that for any L > 
and ^l^l) < E < 1 there exist unique radii < r_(E,L) <r L < r + (E,L) < oo 
such that 

V L (r±(E,L)) = E, and ¥ £ (r) <E <^ r-(E,L) <r<r+(E,L). 

For later use we note that by the implicit function theorem the mapping 
]0,oo[9Li— >tl is continuously differentiable, since 4-(mo(r) — L/\/L + r 2 ) > 
0. The same is true for the mapping (E,L) *-^r±(E,L) on the set {(E,L) e 
Ex]0,oo[| ^ L (r L ) <E<1}, since \l/ L (r)^0 for r^r L . We also note that 

v^(r L ) = 4vrp (r L ) + -j > 0. (3. 16) 

If ft. is to solve Eqn. (J3.13)) . then for any solution r i— > (r(r),u;(r),L) of the 
characteristic system (|3.14j) . 

—h(r(T),w(T),L)=g(r(T),w(r),L). 
As long as w we can rewrite this in terms of the variable r: 



— h(r,w{r),L) = — g[r,w{r),L). 

ar w{r) 
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Since E and L are constant along characteristics we express w in terms of 
these conserved quantities via the relation (J3.15j) . More precisely, with 

E-U (r) 



w(r,E,L) := y/(E - U Q {r)f - L/r* - 1, q(r,E,L) := 01 (3.17) 

w{r,E,L) 

we have w(r) = signw(r) w(r,E,L) and 

—h(r,w(r),L) = sign it; (r)q(r, E , L) g(r,w (r, , L) , L) . 

Outside supp/o w e set /i = 0. For (r,w,L) G supp/ with L>0 we define as 
follows. We let E: = E(r,w,L) so that r_(E,L) <r <r + (E,L), and 

pr 

h(r,w,L) := signw / g(s,w(s,E,L),L)q(s,E,L)ds. (3.18) 

Jr-(E,L) 

In order for this definition to be consistent we show that the integral above 
vanishes if r = r + (E,L) so that h(r,0,L) = 0. To this end let G be as specified 
in Eqn. (|3.6|) so that JJ dfG((f)(E),L)g = 0. We want to use the change of 
variables (x,v) i— > (r,w,L) t— > (r,E,L). Since ^ is even in f and hence in u> we 
can extend this integral only over {w>0}, and 

cfo = 8vr 2 rfr dL = 8n 2 q(r, E,L)drdEdL. (3.19) 

Thus 

»r+{E,L) 




g(r,w(r,E } L) } L)q(r,E,L)drd f G((j)(E),L)dEdL = 0, 

MJr-{E,L) 

where M :— {(E,L)(x,v) \ fo(x,v) > 0}. The class of test functions 
dfG(<f)(E),L) is large enough to conclude that for almost all E and L the 
integral with respect to r vanishes as desired. 

As defined above the function h is not sufficiently regular to apply 
Lemma 13.31 to it. Fubini's Theorem shows that h is measurable, because 
by the change of variables formula the function 

(s,r,E,L) h-> l[y L (r L ),E )(E)l[ r _(E,L),r + (E,L)](s)l[0,r)(s) 

g(s,w(s,E,L),L)q(s,E,L) 

is integrable; r <max{|x| | (x,v) G supp/o} ; and 1m denotes the indicator 
function of the set M. 
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The cut-off h. 

We need to regularize h, and as a first step we study a cut-off version of h. 
To do so we need an auxiliary result. 

Lemma 3.4 For every mGN there exists a constant C m >0 such that for 
L>l/m and ^ L (r L )<E<E , 



rr+(E,L) 

/ q{r,E,L)dr<C m . 

Jr-(E,L) 



Proof. We first claim that for all m e N there exists a constant r\ m > such 
that for allL>l/m, ^ L (r L )<E<E , aadre[r-(E,L),r+(E,L)], 

\%(r)\ 

' ! >r]m- (3.20) 



V*L(r)-^L(r L 
To see this, let mGN and L,E,r be as specified. Then 



#o > E> * r (r) = C/ (r) + V 1 + L /^ 2 > ^o(O) + y/1 + l/(mr 2 ), 

and hence r > (m((E — f/ (0)) 2 — 1)) . Let R := max{|x| | (x,v) G supp/o}- 
Then 

L < r 2 ((£ - U (r)) 2 - l) < i? 2 ((E - U (0)) 2 - l) . 

If the above claim were false, there would exist a sequence (r n ,L n ) — > (f,L) 

in the set [(m((E - U (0)) 2 - 1))~ 1/2 x [l/m,i? 2 ((E - f/ (0)) 2 - 1)], such 
that 

lim - L " V " 7 =0. (3.21) 

We can assume that f—ri since otherwise ^^(f) > ^Z^(r^) and \l/^(f) = 
which is a contradiction to the uniqueness of the minimizer ri of By 
Taylor expansion at r = we find intermediate values 9 n and r n between 
r n and r^ n such that 

\K n (rn)\ M n (0n)(rn-r Ln )\ 



so that 1)3.21)1 contradicts ()3.16)1 . and ()3.20)) is established. 



2\m{r L )\, n-^oo, 
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We split the integral under investigation into two parts, 

r+(E,L) rr L fr+(E,L) 

q{r,E,L)dr = •••+ / ...= : / 1 + / 2 . 

r-(E,L) Jr^(E,L) Jr L 

Using the definitions of w(r,E,L) and and the fact that E — Uo(r)> 
ty L (r)-U (r)>l we find that 

q{ r,E,L) E ~ U ^ E - U ^ 



w(r,E,L) ^/(E-U (r))*-(* L (r)-Uo(r))< 

E-U (r) 



y/{E - V L (r))(E + L (r) - 2U (r)) 



< E-U (0) < C 



'2(£-¥ L (r)) ^-^(r) 

for some constant C which does not depend on r,E, or L. In the inte- 
gral Ji we change variables w = a/^lM — ^l(^l) so that ^f: = ^^lM < 
on [r_(£ , ,L),ri[. By ()3.20j) and the previous estimate, 

h < C f° 1 ~du 

Jy/E-v L {r L ) y/E--$ L {r L ) -u 2 du 

1 rV E -*^) du 1 f 1 ds 

< — / = — / <oo. 

VmJo ^/E-y L (r L )-u 2 VmJo Vl-s 2 

We argue similarly for I2, and the proof of the lemma is complete. □ 
For m e N let 

Q m := \ (x,v) e M 6 I E(x,v) <E -—, L(x,v) > — \ . 
{ mm) 

With the help of the previous lemma we can prove that ln m h is square 
integrable and solves the equation {— E,l^ m h} = lo m <? in the sense of distri- 
butions. 

Lemma 3.5 Form&N large, ln m h 6 L 2 (M 6 ) ; and for any spherically sym- 
metric test function ip — ip(r,w,L) G C 1 ([0,oo[x]R x [0,oo[), 

J J {-E,il)}lu m hdvdx = -JJ ijjl Um gdvdx. 
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Proof. In order to prove that ln m h G L 2 (R 6 ) we observe that the integrand 
is even in v so that we can apply the change of variables (|3.19l) : 



rr rr rr+{E,L) 

l nm h 2 dvdx = 8n 2 / h 2 {r,w{r,E,L),L)q{r,E,L)drdEdL, 

J J J J S m Jr_(E,L) 



where 

S m :={(E,L) = (E,L)(x,v) \ (x,v) G fi m } . 

Let (E,L) G S m . Writing w(r) — w(r,E,L), q(r) = q(r,E,L) and r± — r±(E,L) 
for brevity the definition (j3.18|) of h, the Cauchy-Schwarz inequality, and 
Lemma [3.41 imply that 



r+ pr + 

h 2 (r,w(r) ,L)q(r) dr - 



g(s,w(s),L)q(s)ds 



q{r)dr 



< I / q(r)dr / g (r,w(r),L)q(r)dr 



r+ 



<C m / g (r,w(r),L)q(r)dr. 

J r_ 

A further integration with respect to E and L and the change of variables 
(r,E,L) I— > (x,v) shows that ||ln m /i||2 is bounded in terms of C m and | 1 1 2 . 

Let ip be a test function as specified in the lemma. As above, we param- 
eterize the solutions of (j3.14j) by r, distinguishing between w>0 and w<0. 
Then it is easy to show that 

{-£,V>} = sigmo-^— [ip(r,w(r),L)). 
Using the same change of variables as before and recalling Eqn. (|3.18J) we 
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find that 



J J {~E,ip}lQ, m hdvdx- 



••• + 

{w>0} J{w<0} 



In 2 II I —^—ty(r,w(r),L)]h(r,w(r),L)q(r)drdEdL 
q(r) dr 



. , , , i r+ 1 d 



Sir // / —— — [ijj(r,—w(r),L)]h(r,—w(r),L)q(r)drdEdL 
JJsm Jr^ <?( r ) dr 

-8n 2 / ip{r,w{r),L) g(r,w(r),L)q(r)drdEdL 

r r r r + 

2 



7T / / / ip(r,—w(r),L) g(r,—w(r),L)q(r)drdEdL 

J J Sm J T- 

4>ln m gdvdx; 

note that h(r,±w(r),L) = for r = r±(E,L) which together with the defini- 
tion of h is used in the integration by parts above, and also that g is even 
both in v and w. The proof of Lemma [3.51 is complete. □ 

Regularization of h. 

The function ln m h is not smooth, so in order to apply Lemma EO we smooth 
it. For mGN fixed the function ln m h function of r,w,L, supported 

in a cube of the form 



Q-^lRoMxl-WotWolxlLotL 



ih 



with 0<.Ro<.Ri, W >0, and 0<L <£i- Let (eC™(R 3 ) be even in 
all three variables, i.e., ((zi,z 2 ,Z3) = ((\zi\, \z 2 \, \zs\), (>0,f( = l, and let 
Cn( z ) :=n 3 ((nz) for For n sufficiently large we define 

/»oo /»oo /*oo 

h n (r,w,L):= / / (l nm h)(f,w,L)( n (r-f,w-w,L-L)dLdwdf. 

JO J-odJo 

Clearly, ft, n G C^°(]0,oo[xIRx]0,oo[), and without loss of generality we can 
assume that supp h n C Q H {E < E } . Since /i is odd in w, so is h n , and 
clearly /i n — > ln m h in L 1 flL 2 . The crucial step is to show that 

\im{-E,h n } = l nm g inL 2 . (3.22) 
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We fix (r,w,L), write p: = ^Jl + w 2 + L/r 2 , p: = \Jl + w 2 + L /r 2 , and / :- 
Jo I-oo Jo ^ or brevity, and split the convolution integral as follows: 



w 



{-E,h n } = —d r h r , 
P 



L 



r 3p 



U' Q {r)\d w hr 



ln m h){r,w,L) 



[ln m h)(r,w,L) 



V V J 



Cn{ r — r,w — w,L — L)drdw dL 



w 
—i 
P 



L 



£ n (r — r,w — w,L — L)dr dw dL 

=:Ji + J 2 - 
According to Lemma f3. 51 

h = -J ln m h{-E,(n(r-;w-;L--)} 

= J ( 1 nrn9)(r,w,L)( n (r-f,w-w,L-L)dfdwdL^ln m g in L 2 . 

We show that Ji converges to as n— >oo. To this end we introduce 
new variables r = n(r — f), w = n(w — w), L = n(L — L). Then f = r — r/n, 
w = w — w/n, L = L — L/n, dfdwdL = n~ 3 dfdwdL, and d r ( n (r — f,w — w,L — 
L) =n 4 df((r,w,L). Now we note that 

(w w\ _ _ _ _ - 

— ) df( n (r — r,w — w,L — L)drdwdL 

The analogous identity holds for the second part of J 1( and we find that 

u> — w 



w w 
n | — 



n- 



P 



1 1 
+ nw | 

J> P. 
w 2 — w 2 + L/r 2 — L/f 2 



w 

\-nw 

p pp{p+p) 

w 2ww 2 + wL/r 2 — 2rwL/r 3 

p 2p 3 
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since 

L L L — L (f—r)(f+r) 

We easily see that 



Furthermore, 



n{U' {T)-U' {f))^U'^r)f. 



L \ { L/f — L/r L (I 1 



r 6 p r 6 p j \ p r 6 \p p 



1 / 3Lr L\ L 2ww + L/r 2 -2rL/r 3 



Thus most of the terms above integrate to in the limit of the expression 
Ji, leaving us in the end with 

Liu f 

■h^(^n m h){r } w } L)-^-p J (-fd f + wd^)C(f,w,L)dfdwdL = 0, 

where the limit is in L 2 , and arises in the last equation from an integration 
by parts with respect to f and w and an exact cancellation. This completes 
the proof of 

The contradiction to Lemma \S. <A 

The functions h n have all the properties required in Lemma By 
Eqn. ()H.22j) . lim n ^ 00 {-E, h n } = l Um g in L 2 , and since h n and ln m g are sup- 
ported in a common compact set, VU{_E,h n } ~^ VC/i„ g in L 2 as n —>■ oo. This 
implies that 

D 2 H c (fo) [{-E,h n }} -> D 2 H c (fo) [W], n^oo. 

Since lim m _^ 00 lQ m (7 = (7 in L 1 flL 2 (IR 6 ), there exists m sufficiently large such 
that ln mo 9 7^ 0, and by Lemma 1331 la TO0 /i^ 0. For all m>m , Q mo C f2 m and 
thus ln mo ^ 2 < ln m h 2 . Since <fi'{E) < and [/^(r) > we have for all m > m , 

^c(/o)[la^]>-^ // ^yi^la mo ^(l + |^ 2 )^ 3/2 ^ C /x>0. 

Letting m — ► 00 we get a contradiction to Lemma 13.21 and the proof of The- 
orem is complete. 
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4 Global existence of the perturbed solutions 



Let fo be a steady state as specified in Section 2, and choose 5 > according 

o 

to Theorem 12 .21 For solutions with initial data fET>f satisfying the estimate 

o 

d(f,fo) < 5 the potential energy and therefore also the kinetic energy remain 
bounded. In the present section we show that this implies that such solutions 
are global in time. 



Proposition 4.1 Let /gD/ . Then there exists a unique, smooth solution 

o 

/ of the relativistic Vlasov-Poisson system with f(0) = f on a maximal time 
interval [0,T[ with T>0. If 



then this solution is global, i.e., T = oo. 

Proof. It is well known that a compactly supported, non-negative initial 

o 

datum /eC*(R 6 ) launches a unique C 1 solution on some time interval [0,T[ 
which we choose maximal. A detailed proof for the non-relativistic case 
is given in 20], and the proof carries over to the relativistic situation. In 
addition, the solution is known to be global if 



However, there is a technical problem here. The steady state fo and therefore 
also data taken from T>f need not be continuously differentiate on IR 6 , they 
need only be continuous on R 6 and continuously differentiable on the set 
where they are strictly positive. But due to the spherical symmetry of the 
functions in Vf these regularity properties propagate to the solutions, since 



is continuously differentiable if / and therefore p are continuous. In par- 
ticular, the characteristic flow of the Vlasov equation is well-defined and 
continuously differentiable and / is constant along these characteristics. To 
turn these arguments into a rigorous proof would require going through the 
local existence proof for the present situation, which would be lengthy but 
straight forward. 



SUp #kin(/(*))<00 



0<t<T 



sup{|f | | (x,v) € supp/(£), <t < T} < oo. 
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Let us now assume a bound on the kinetic energy of a maximal, spherically 
symmetric solution / on some time interval [0,T[. We want to deduce a bound 
on the quantity 

P(t) := sup{|i;| | (x,v) G supp/(s), < s < t}, 

which by the continuation criterion recalled above implies global existence. 
The arguments which follow are adapted from [3] , but we prefer to give a self- 
contained presentation. First we note that a standard interpolation estimate 
implies that ||p(t)||4/3 is bounded on [0,T[, cf. j2Dl Lemma 5.1]. This in turn 
implies that for all < s < t < T, 

1 1 VC/( S ) | U < C| 1 1^ 9 < CJ- 5 / 3 (t), 

o 

cf. [213 Lemma PI]; the constant C depends only on the initial datum / and 
on the bound on the kinetic energy. Due to the above formula for the field 
in the spherically symmetric case, 

\U'{s,r)\ <cmin{r- 2 ,P 5/3 {t)}<c(r + a)- 2 , r>0, (4.1) 

where we abbreviate a := P~ 5//6 (t); the constant c>0 is now kept fixed. 
For L > we define the function 

£ L {u):=Vl+u 2 -^ , u>0. 

yL + au 

It is easy to check that £l g C 2 ([0,oo[) with £^>0. Moreover, there exists a 
unique ul > such that ^l(ml) = 1 and £l is strictly increasing on the interval 
[ul, oo[. Clearly, 

Zl(u)>Uu) = VT+^-- 

a 

for all u > and L > 0, and this implies that 

u L <uo = V (1 + c/a) 2 - 1. 

Consider now a characteristic [0,T[3 s i— > (x(s),v (s)) in the support of /. 
With 

r(s) :=\x(s)\, u(s) := \v(s)\, w(s) :=x(s) -v(s)/r(s), L(s) := \x(s) x v(s)\ 2 
it follows that L(s) = L is constant and 

^-u 2 (s) = -2U'(s,r(s))w(s). 
as 
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Let us first consider the case that L>0, fix some time i£]0,T[, and assume 
that w(t) < 0. We choose < t < t minimal with the property that w(s) < 
for s&]to,t]. Then u(t)>u(t ), and 

= (tt^-u'Ms))) /°° <o 

\(r(s)+a) 2 "J y/l+ u \s) 

by the estimate (j4. ljl . Since u 2 (t) =w 2 (t) + L/r 2 (t) we have the estimate 
r(i) > vL/u(t), and hence 

, r°° dr , r°° dr 

Jr(to){ r + a ) JVL/u(t)(r + ay 

If w(to) =0 then u(t ) = vL/r(t ) and hence 

&(«(*))<&(«(*))), 

which by choice of it^ and the fact that >«(to) implies that -u(t) <ul< 
uq. If instead w(to) < then to = and < a/1 + P 2 (0). Hence in the 

case L > and w (t) < we have shown that 

^/l+u\t)< v / l + P 2 (0) + c/a; (4.2) 

notice the formula for ito stated above. 

Next we consider the case that L>0 and w(t)>0. Again, we choose 
< to <t minimal such that w(s) > for s 6 [to, t\. Then u(t) <u(t ). If t = 
then u(t)<P(0). If t >0 there must exist a sequence 0<tj/*t sucn that 
w(ti) <0, which by the case we already dealt with implies that 

^l+u 2 (U)< v /l + P 2 {0) + c/a, 

and by continuity, this estimate must also hold with to instead of £j. Hence 
the estimate ()4.2j) holds for all characteristics with L > 0, but since the initial 
data of a characteristic in supp/ with L = can be approximated as closely 
as we wish by data with L > and since characteristics depend continuously 
on their initial data the estimate (|4.2|) remains true if L = 0. Hence 

a/1 + P 2 (t) < a/1 + P 2 (0) + c/a = a/1 + P 2 (0) + cP 5/6 (t) , 0<t<T. 
This implies that P is bounded on [0,T[ by some constant which depends on 

o 

the initial datum / and the bound for the kinetic energy, and the proof is 
complete. □ 
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5 Stationary solutions 



As we explained at the beginning of Section 2 steady states of the relativis- 
tic Vlasov-Poisson system can be obtained by making the ansatz (|2.1|) for 
some function <p — (p(E). This ansatz solves the stationary Vlasov equation 
with potential Uq in the sense that fo is constant along characteristics. By 
substituting it into the definition of the spatial density the latter becomes a 
functional of the potential Uq: 

po(x)=g < f > (U (x)), (5.1) 

where 

poo 

^(w):=4vr / (p(E)(E — u) \{E — u) 2 — l] dE, uGl. (5.2) 

Ju+l 

Hence a steady state of the relativistic Vlasov-Poisson system is obtained, if 
for a given choice of <fi it can be shown that the semilinear Poisson equation 

A[/ = 47T^([/o) 

has a solution on M? such that f$\=(f)oE has finite mass and compact sup- 
port. Since in this section we consider only stationary solutions, we drop the 
subscript and denote the steady state by / and its potential by U . 

One can show that any solution obtained by the above approach must 
a-posteriori be spherically symmetric, so we only look for solutions U = 
U(r), r >0, of the equation 

An r 

U'(r) = — j 9(i> (U(s))s 2 ds,r>0. (5.3) 

In |2j, Batt showed that the polytropic ansatz 

<f>(E) = (E -E) k + , 

with fce]0,2[ leads to steady states with compact support and finite mass. 
Below we present two approaches which can be used to prove the existence 
of steady states for more general functions 0. First however, we show that 
4> must vanish for large values of E if the resulting steady state is to have 
finite mass. 
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Proposition 5.1 Let >[0,oo[ be measurable and let (f,U) be a sta- 

tionary solution of the relativistic Vlasov-Poisson system in the sense that 
f(x,v)—<p(E) and U G C 1 ([0,oo[) solves \5.3\) with given by \5.ty) . Let 
JJfdvdx < oo. Then C/qo := lim^oo U{r) < oo and (fi(E) = a.e. for E> 

Uoo + 1. 

Proof. Let M : = jj fdvdx. Since by (Ql, < U'{r) < M/r 2 , U is increasing 
with a finite limit Uoo for r^oo, and since is decreasing, 

poo poo 

M = 4n g<f> (U(r))r 2 dr>4n g^U^r 2 dr, 
Jo Jo 

which means that g^ ) {U OQ ) = Q. This is only true if <p(E) = for almost all 
E>U 00 + 1. □ 

In what follows we are only interested in steady states with finite mass so we 
always assume such a cut-off energy E . 

Under suitable assumptions on <ft Eqn. (J5.3j) has a unique solution on 
[0,oo [, if we prescribe U(Q). 

Proposition 5.2 Let 0:IR^ [0,oo[ have the following properties: <fi is mea- 
surable, there exists a cut-off energy E ^M. such that (f>(E) = for E>E , 
and there exists k>—l such that for every E\<Eq there exists a constant 
C > such that 

(f>(E)<C(E -E) k , Ee[E u Eo[. 

Then for every u EM. there exists a unique solution U GC 1 ([0,oo[) of \5.3\) 
with U(0) =u . 

Proof. One can show that under the given assumptions on <ft the function 
g^ defined in Eqn. (j5.2|) is continuously differentiable. The existence and 
uniqueness of U on some interval [0,5] then follows by a contraction argument. 
Since U is increasing, either U < Eq — 1 on its maximal existence interval, and 
the solution is global, or there exists some R> such that U(r) > Eq — 1 for 
r>R, which implies that g^(U(r)) = for r>R, and again U is global. □ 

Remark. The solution obtained above will in general not satisfy the bound- 
ary condition = which we required in previous sections, and which is 
part of our formulation of the relativistic Vlasov-Poisson system. However, 
the solution has a finite limit 1/^, and by subtracting this limit from U and 
redefining the cut-off energy Eq accordingly we obtain a steady state with 
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the right boundary condition. It is clear that E and U(0) cannot be freely 
prescribed under the given boundary condition for U at infinity. 

We now proceed to our first result which guarantees the existence of 
compactly supported steady states with finite mass. We follow an approach 
which was used for the Vlasov-Poisson and Vlasov-Einstein systems in [22] • 

Proposition 5.3 Let <fr be as in Proposition \b."A and assume in addition 
that 

<p(E) = c(E -E) k + O((E -E) k+s ) as E^E -, 

for some — l/2<k<3/2, c>0, and 5>0. Let (f,U) be an induced steady 
state of the relativistic Vlasov-Poisson system, i.e., f(x,v)=<f)(E), and U G 
C 1 ([0,oo[) satisfies Then this steady state has compact support and 

finite mass. 

It is remarkable that up to technical assumptions only the behavior of at 
the cut-off energy E needs to be specified. In particular, there are plenty of 
functions <ft which satisfy the assumptions of this theorem and of the stability 
result in Section 2 as well. Beside the polytropic states with l<A;<3/2 a 
notable example which frequently appears in the astrophysics literature is a 
King type model 

f(x,v) = <P(E) = (e E °- E -l) + . 

The basic set-up of the proof of Proposition 15.31 is as follows. If U (0) > 
Eq — 1 then this relation holds also for r > 0, and the steady state is trivial. 
Hence we consider a solution U such that U(0) <E — 1, and we define [0,R[ 
as the maximal interval on which U < E — 1 so that supp p=[0,R}. If R = oo 
then Uoo = lim r _ >00 L r (r) < Eq — 1 exists by monotonicity. Now assume that 
Uoo <Eq — 1. Then p(r) > ^g^iUoo) > which implies that U'(r) > Cr with 
C > which upon integration gives the contradiction that Uoo = oo. Hence 

U(r)<E -l on [0,R[, and lim U(r) = E -l; (5.4) 

if R < oo the limit assertion is true as well. In order to show that (J5.4j) indeed 
implies that R < oo we use a special case of a result due to M akino [13] ; a 
proof can also be found in |22l Lemma 3.2]. 

Lemma 5.4 Let x,y G C 1 (]0,-R[) be such thatx,y>0 and 

rx' = a(r)y — x + x 2 
r v' — V ( a — /3(r)x) 
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on]0,R[, where a>0 is a constant, a, [3 G C(]0,R[), inf r . g ] 0R [Q;(r) > 0, and 
\im r ^ R _P(r) e]0,a[. Then R< oo. 

Proof of Proposition 15.31 We first introduce the radial pressure 

/fx ' v \ ^ 
y j f(x,v)dv. 

Like the spatial density p the pressure becomes a functional of U via 
Eqn. 

p(r) = ^(t^(r)) 

where 

h^ u ):=— [ <j){E){E-u) [(E-u) 2 -l] 3/2 dE, ueR. (5.5) 
3 Ju+i 

On the interval }0,R[ we define 

rj(r) :=E — U{r) — 1, m(r) : = 47r / s 2 p{s)ds 

Jo 

and 

z(r ):=— y(r):=Aixr 2 ^- . 
rr]{r) p[r) 

A short computation shows that the latter functions satisfy a system of ODEs 
of the form stated in Lemma f5.4l with a:=2 and 

p h(jj) 



a 



pv g{v)v 



p p J m g(rj) h(rj) 
For the sake of convenience we have redefined 

g{ V ):= 4vr [\( e )(l + 7i- e) [2{ V -e) + {v~ rf] ^ de, 
Jo 

4-TT P q /n 

h{rj):=-J^ ^e)(l + ^-e)[2^-e) + (^-e) 2 ] 3/2 de J 

and 

^(e): = 0(E o -e) = ce fc + O(e fc+5 ) as e = E -E^0 + . 
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In order to investigate the asymptotic behavior of (5 we also need 

</fa)=47T rV( e )[(2(^- e ) + (r ? - e ) 2 )- 1/2 + 2(2(r ? - e ) + (r ? - e ) 2 ) 1/2 l de, 
Jo L J 

/ i '(r ? )=47r^%(6)[(2(r ? - e ) + (r ? - e ) 2 ) 1/2 + l(2(r ? - e ) + (^- e ) 2 ) 3/2 



cie. 



We need to substitute the asymptotic expansion of ip into g, h, and their 
derivatives; note that < rj(r) — > as r — > i?— . The leading order terms always 
contain an integral of the form 



f e k {r 1 -e) m de = r 1 k+m+1 [ s k (l-s) m ds=:r] k+m+1 c k ^ 
Jo Jo 



m> —1, 



where the constants c k ^ m satisfy the relation 

c k , m ~i _ k + m+1 



Hence for 77 — > 0+, 



, m>0. 



g(r]) = 47rc2 1 / 2 c kA/2 r ] k+3 / 2 + 0(r ] k+ V 2+s ), 

M^) = yc2 3 / 2 c fc ,3/ 2 ^ +5 / 2 + 0(^ 5 / 2+5 ), 

g\r ] )=Anc2-^c k ^ /2 r ] k+ ^ + 0{r ] k ^), 
h'{rj) = Anc2 1 / 2 c kA/2 r ] k ^/ 2 + 0(r ] k ^ 2+s ). 

If we substitute these expansions into the formula for a we find 
a{r) = l^/l + 0( V s (r)) - (k + 3/2 + 1)' 1 > 0, r - R-, 

O C k ,l/2 

Clearly, a(r) > for < r < R and since 77(0) > 0, 

lim a(r) = , r) J / > 
1 } 77(0)^(77(0)) 

so that inf J . 6 ] 0iR [Q;(r) >0 as desired. As to /3, 
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and by the assumption on k this limit lies in the interval ]0,2[. We can 
therefore apply Lemma f5. 41 to conclude that R< oo, and the proof of Propo- 
sition is complete. □ 

Although the class of steady states obtained by the approach above is 
quite large, the condition on k excludes for example polytropes with k>3/2 
and possible generalizations of those. Hence we present a second approach 
the idea of which is the following. If one introduces the speed of light c as a 
parameter into the relativistic Vlasov-Poisson system then solutions of this 
system tend to solutions of the non-relativistic case as oo. Exploiting this 
for Eqn. (|5.3|) one can show that the potential in the relativistic case has to 
cross the cut-off energy for sufficiently large c, if the corresponding potential 
in the non-relativistic case does. In this way one can show that essentially 
every compactly supported steady state of the non-relativistic Vlasov-Poisson 
system has its relativistic counterpart. In particular, this will be true for the 
polytropic ansatz with 0<k<7/2, which extends the range found in [2J. 
The perturbation argument outlined above was used for the Vlasov-Einstein 
system in |21j . 

If the speed of light is not normalized to unity in the relativistic Vlasov- 
Poisson system the only change is that a/1 + \v | 2 must be replaced by 
\A + \v\ 2 /c 2 . In the sequel we use the parameter 7:=l/c 2 . In order to 
obtain the correct non-relativistic limit as 7— > we adjust the ansatz (|2.1j) : 

f(x,v)=<p(-y/l + 'y\v\ 2 + U(x) - - J . (5.6) 
V7 7/ 

This ansatz satisfies the Vlasov equation with potential U and parameter 7. 
As before the spatial density p becomes a functional of the potential U, 

p(r)=94,(j,U(r)), 

where 

POO 

^(7,m):=4tt / <j>(E) (l + 7 (£- M )) [2(E -u) + 1 (E -u) 2 ] 1/2 dE, uER, 

J u 

and it remains to analyze the Poisson equation for U which takes the form 

47r r 

U'(r) = — g^,U{s))s 2 ds,r>U. (5.7) 
r 2 Jo 
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For 7 = this is precisely the equation which arises if one looks for steady 
states of the non-relativistic Vlasov-Poisson system via the ansatz 

f(x,v) = (f>(~\v\ 2 + U(x)\ (5.8) 

and an analysis of the limit 7— >0 yields the following result. 

Proposition 5.5 Let <f) E L™ C {R) , <f)>0, and <f>{E) = 0, E>E . Let (f ,U ) 
be a non-trivial, compactly supported steady state of the Vlasov-Poisson sys- 
tem obtained by the ansatz \5.8\) . i.e., Uo E C 1 ([0,cxo[) solves Eqn. \5. l\j with 
7 = 0, and Uq(0) < E but Uq(R) > E for some R>0. Then for all 7 > suf- 
ficiently small the solution U EC 1 ([0,oo[) of \5. l\j with £7(0) — Uq(0) satisfies 
U(R)>Eq, i.e., the resulting steady state of the relativistic Vlasov-Poisson 
system with parameter 7 has compact support and finite mass. 

Proof. It is easy to see that there exists a constant C > such that 

\g l j > ('y,u)-g(0,uo)\<C'y+C\u-uo\, u,u >U (0). 

In what follows, constants denoted by C are positive, may change from line 
to line, and never depend on 7, U, or r. By monotonicity the solutions U 
and Uq under consideration take only values u>U (0), and hence 

\U'(r)-Uo(r)\<Cj+C / \U(s)-U (s)\ds, 0<r<R. 

Jo 

If we integrate this estimate, observe that U (0) = U o (0), and apply Gronwall's 
lemma we find that 

\U(r)-U (r)\<C^, 0<r<R. 

The latter holds in particular for r — R, and since by assumption Uo(R) > E Q 
the same is true for U, provided 7 is small enough. □ 

Remark. The requirement that 7 must be sufficiently small in the re- 
sult above actually changes the relativistic Vlasov-Poisson system we started 
with. However, it is easy to check that if / is a steady state of the system 
with some 7>0 then the rescaled function f^(x,v) :=7~ 3/ ' 2 /(7~ 1//2 x,7~ 1 / 2 f ) 
is a steady state of the original system with 7 = 1. This scaling changes 
the cut-off energy, which however has to be adjusted in order to sat- 
isfy the boundary condition U(oo) = anyway, and it changes the ansatz, 
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/ 7 (a;,f)=7 3,/2 0((a/1 + \ v | 2 + U(x) — 1)/7), but it does not change the prin- 
cipal way in which the ansatz depends on the particle energy. 

In the non-relativistic case it is well known that the polytropic ansatz 
with 0<k<7/2 leads to compactly supported steady states; as a matter of 
fact the compact support is lost for k — 7/2, the so-called Plummer sphere, 
and finite mass is lost for k>7/2. Generalizations of these polytropes were 
obtained by variational techniques in [TU] , and by techniques which are closer 
to the ones used to prove Proposition 15. 31 in so there are many examples 
which one can extend to the relativistic case by Proposition 15.51 

To conclude this section we emphasize that its purpose was only to demon- 
strate that there are many examples of steady states which satisfy the as- 
sumptions of our stability analysis. Its purpose was not a systematic investi- 
gation of the possible steady states of the relativistic Vlasov-Poisson system. 
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